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ABSTRACT 



We present Wide Field Camera 3 images taken with the Hubble Space Tele- 
scope within a single field in the southern grand design star-forming galaxy M83. 
Based on their size, morphology and photometry in continuum-subtracted Ha, 
[S II ], H/3, [O III ] and [O II ] filters, we have identified 60 supernova remnant 
candidates, as well as a handful of young ejecta-dominated candidates. A cata- 
log of these remnants, their sizes and, where possible their Hct fiuxes are given. 
Radiative ages and pre-shock densities are derived from those SNR which have 
good photometry. The ages lie in the range 2.62 < log^Tj-ad/jT^) < 5.0, and the 
pre-shock densities at the blast wave range over 0.56 < ng/cm^^ < 1680. Two 
populations of SNR have been discovered. These divide into a nuclear and spiral 
arm group and an inter-arm population. We infer an arm to inter-arm density 
contrast of 4. The surface fiux in diffuse X-rays is correlated with the inferred 
pre-shock density, indicating that the warm interstellar medium is pressurised by 
the hot X-ray plasma. We also find that the interstellar medium in the nuclear 
region of M83 is characterized by a very high porosity and pressure and infer a 
SNR rate of one per 70-150 yr for the nuclear {R < 300 pc) region. On the basis 
of the number of SNR detected and their radiative ages, we infer that the lower 
mass of Type II SNe in M83 is M^in = 16^1 M©. Finally we give evidence for 
the likely detection of the remnant of the historical supernova, SN1968L. 

Subject headings: supernovae: general-ISM: structure, supernova remnants- 
galaxies: ISM, starburst, structure 



1. Introduction 



The techniques for the detection of supernova remnants (SNR) in external galaxies were 



pioneered by Mathewson & Healey (1964); Westerlund & Mathewson (1966) and Mathew 



son & Clarke (1973). These techniques used the characteristic non-thermal radio spectrum 



produced by SNR along with their optical characteristics, especially the relative strength 



of their [S II ]AA6717, 31 lines relative to Ha. Later, Dopita (1977a) established that these 
optical characteristics were a natural consequence of the radiative shock waves associated 
with the propagation of the supernova blast wave through the interstellar medium (ISM) 
of the host galaxy, and that the details of the radiative spectrum in turn depends on the 
chemical abundances in the ISM of the host galaxy, rather than on the chemical abundances 



in the stellar ejecta (Dopita 1977b). Strong [S II ]AA6717, 31 and [O II ]AA3727, 9 emission 
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arises in the recombination zones of these radiative shocks, where much of the Balmer hne 
emission is also produced. 

The characteristic strength of the [S II ]/Hq; ratio seen in SNR (typically greater than 
0.4) has proven to be an extraordinarily versatile technique for detecting regions excited by 
radiative shocks. The association of such regions with either non-thermal radio emission or 
X-rays emitted from the shock-heated gas behind the SNR blast wave has usually sufficed to 
confirm the identification of a SNR. This technique has been successfully used to survey for 



other local group galaxies ( 


D'Odorico, Dopita & Benvenuti 


1980 


), and in M33 ( 


D'Odorico 


&; Sabbadin 


1978 


Blair, Kirshner & Chevalier 1981 


Long et al. 


1990 


Gordon et al. 


1998 



Long et al. 2009) 



M83 (NGC 5236) is a southern Grand-Design SAB(s)c galaxy. It has an inclination angle 
of 24 degrees and has an angular size of ~ 13 '. Its distance is given in the NASA/IPAC Ex- 
tragalactic Database (NED) as 4.56 ±0.09 Mpc using a combination of redshift-independent 
techniques. It has well-formed spiral arms and distinct dust lanes. Detailed observations at 
radio ( |Crosthwaite et al.||2002[ [Maddox et aL]|2006D , infrared ( fVogler et al.|[2005| [Rubin et al." 
2007| , ultraviolet dBohlin et al.|[l990j [[Boisser et al.|[2005| [Thilker et al.||2005D , X-ray ( [Soria fc 
Wu||2002| |2003| |Kilgard|2005D and optical wavelengths ( [Calzetti et al.|2004D indicate that the 
nucleus and spiral arms of M83 are regions of intense star formation. [DOdorico, Pettini, &: 



Ponz ( 1985 ) found absorption features in the nucleus consistent with a population dominated 



by massive OB stars, which are believed to be the precursors of Type Ib,c/ II SNe (Huang 



1987). Observations by Jensen, Talbot & Dufour (1981) revealed a large number of young 



massive stars along the inner arms, indicating active star formation in these regions as well. 
These observations are consistent with a recent and ongoing burst of star formation lasting 
over at least the last ~ 10'^ (Trinchieri, Fabbiano, & Palumbo 1985). Indeed, it has been 



the site of no less than six optically-detected supernova events in the past century making it 
the most active galaxy in this respect along with NGC 6946. This ubiquitous star formation 
and its high surface density makes M83 a natural target for SNR studies. Indeed it has 



been the subject of a detailed study by Blair &; Long (2004), who identified a total of 71 



candidates over the whole face of M83 using the ratio of the [S II ]AA6717, 3lA fines to Ha, 
and followed up with spectroscopy of 25 of these. They found that 15 of their candidates 



were also associated with soft X-ray sources in the catalog of Soria & Wu (2003) or with 



non-thermal radio sources ( Maddox et al. 2006 ) . 



In this paper, we present new observations of a single field in M83 taken by Wide 
Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). These images were obtained 
in many narrow band and broadband filters, and we use these data to discover new SNR 
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candidates not only using the traditional [S II j/Ha ratio, but also by the relative strength 
of their [O II ]AA3727, 9A emission. Twelve of the previously discovered Blair & Long (2004) 
SNR candidates lie within this field. The complete list of candidates is presented in [3| and 
in section |4] we present the interpretation of these data on the basis of their measured sizes 
and Ha fluxes. At the distance of M83, the pixel size of WFC3 (0.0396 arc sec) corresponds 
to a spatial scale of 0.9 pc, and the spatial resolution at Ha (0.06 arc sec) to a distance of 
1.7 pc. This ensures that even the smallest of our SNR candidates is adequately resolved 
spatially. 



2. WFC3 Observations 

2.1. The WFC3 dataset 

The WFC3 data were obtained as part of WFC Science Oversight Committee (SOC) 
Early Release Science (ERS) program (program ID11360: Principal Investigator: Robert 
O'Connell, The University of Virginia). The M83 field, centered at RA=13:37:04.42 Dec=- 
29:51:28.0 (J2000) was chosen to cover the nuclear region as well as the inner spiral arm 
which extends out from the northern side of the nucleus and curves away towards the north 
east. A second, outer field will be obtained later as part of this program. Either three or 
four dithered exposures were taken in each narrow-band filter to remove cosmetic defects 
and assist in the cosmic ray removal. For this paper a sub-set of the data taken in the UVIS 
channel has been used. The observing log for these observations is given in Table [Tj 



2.2. Reduction of Observations 



Following acquisition, the data was passed through the MultiDrizzle software (Fruchter 



fc Hook|[2009| to produce a cosmic-ray cleaned, combined image using the "drizzle" software 



originally described by Fruchter & Hook (2002 ). This delivers a combined, cosmic ray cleaned 



and distortion corrected data product with a built-in world coordinate system. 

In these images, there is a small offset between images taken at different visits to the 
same target. These were removed using the imalign task in the IRAF package. Two weighted 
wide-band images were combined and scaled to provide a continuum reference image for each 
narrow-band image; F336W and F438W for the [O II ] data, F438W and F555W for the B.p 
and [O III ] images and F555W and F814W for the Ha and [S II ] images. The weighting 
of the broadband images was determined by the offset of their bandpass in logarithmic 
wavelength with respect to the emission line in question. It should be noted that the F555W 
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filter also transmits the [O III ] emission line. However, in M83 the strength of this hne is 
generally very feeble. This results from the super-solar chemical abundances in the inner 
region of this galaxy, and the resulting low electron temperatures in the H II regions. As a 
consequence the photometric effect of this leakage is very small. 

The mean redshift of M83 is only 513 km s~^. This is sufficiently small to ensure that 
the line emission is not significantly shifted from the peak of the filter bandpasses. The 
narrow band filter bandpasses are also sufficiently wide to ensure that the likely velocity 
dispersion of the SNR shell (~ 200 — 400 km s~^) is also comfortably accommodated by the 
filter. 

Each scaled continuum image was subtracted from the relevant line image to produce 
a pure emission line image. The subtraction is not perfect, because the color terms for 
the relatively unreddened hot blue stars in the spiral arms are quite different from those 
of the cool, heavily reddened, old stellar population in the nuclear region. The scaling of 
the continuum images was adjusted to provide the best compromise. Under-subtracted or 
over-subtracted residuals are still present, however. 

3. SNR Candidates 
3.1. Disk SNR 

In what follows, we will distinguish between the supernova remnant (SNR) candidates 
in the disk region of M83, and the nuclear region characterized by very high star formation 
rates and generally large reddening. To provide a physical definition of what constitutes 
the nuclear region, we draw a circle with a radius of 300 pc centered at the highly-reddened 
nuclear cluster which we find at coordinates RA=13:37:00.871 Dec= -29:51:55.97 (J2000). 

The selection techniques to find SNR candidates are somewhat different between the 
disk and the nuclear region. Candidates in the general field were identified by blinking 
the aligned and match-scaled Ha, [S II ] and [O II ] images. The SNR candidates appear 
relatively bright on both the [S II ] and [O II ] images compared with Ha. For these 
candidates, we could check the SNR identification by direct measurement of the ratio of [S II 
] to Ha taking into account the fact that the HST F657N filter also transmits the [N II ] 
A6584 hne. The "classical" definition of a SNR is [S II ]/Ha > 0.4. Since the [N II ] hne 
may be as strong as ~ 50% of Ha, we relax the above criterion for identification of a SNR 
to the condition Ff673n/-^F657n > 0.3. Because of their very red colors, old clusters may 
sometimes masquerade as SNR candidates. However, these may be simply removed on two 
considerations - first, they are centrally concentrated, while the SNR candidates appear as 
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shells and second, unlike bona fide SNR candidates, these clusters are very faint in the [O II 
] filter. 

Following identification, the diameters and the positions were measured. The accuracy 
of the coordinates was checked by comparing the positions of the twelve candidates which 



are in common with the catalog of Blair & Long (2004). The agreement is excellent. 



The resultant catalog of the 40 SNR candidates found in the disk is given in Table [2j 
The objects are named by galaxy, nature and field {viz. M83-SNR-1-) and are listed in order 
of increasing RA. The more interesting characteristics of the candidates are given in the 



notes. Several of the candidates have already been identified as X-ray sources (Soria & Wu 



2002| [2003] |Kilgard|[2005| or radio sources ( |Maddox et al.|2006D . Only one ( SN1968L) of the 



six historical SNe listed in Blair & Long (2004) lies within our field. This case is discussed 
in more detail below. 

The Ha images of the disk SNR candidates are shown in Fig [TJ Here, the images 
are shown on a square root stretch to bring out the fainter features. In all cases these 
SNR candidates are clearly resolved, and details of their internal morphologies are evident. 
The sizes of the SNRs are rather similar because all of these candidates are seen in their 
radiative phase. To first order their expansion velocities should be quite similar - governed 
by the timescale for the shock ISM gas to produce a fully radiate shock, and the size will be 
moderated only by the density in the surrounding ISM and the explosion energy. 

Note that, in our images, we find several H II regions with large central cavities. None 
of these have a detected SNR lying in them. This supports the idea that the collective effects 
of stellar winds and ablation of un-ionized inclusions by photoionization has effectively swept 
these cavities clear of H I so that radiative SNR cannot be observed because the cooling 
timescale of the hot, tenuous ISM is too long. 

A few of these SNR candidates are worthy of special mention. Especially interesting is 



M83-SNR-1-14 alias BL39, Blair & Long (2004). This has a linear bipolar structure with 



two bright knots at either end of the long axis, and a more diffuse nebulosity in between, see 
Fig[Tj Its structure is similar to a double-lobe radio source, and is reminiscent of the famous 



galactic supernova remnant W50 (Zealey, Dopita & Malin 1980) associated with the black 



hole binary SS433 (Margon 1979a; Margon et al. 1979b Blundell & Bowler 2004 2005). We 



may speculate that M83-SNR-1-14 too has a black hole binary at its core. BL 39 was one 



of the objects for which Blair & Long (2004) obtained spectra, but apart from the relatively 



weak [O III ] emission, it is not remarkable in its spectral characteristics. Higher resolution 
spectra would be useful to establish whether there is any evidence for a bipolar outflow in 
this object. 
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One might expect that an SS433-hke source should be detected as an X-ray source. 



However, M83-SNR-1-14 (BL39) was not detected with Chandra in X-rays by Soria & Wu 
(2002, 2003) down to the brightness hmit of their survey; 3 x 10^^ erg s~^. However, Namiki 
et aL (2003) give the "typical" 1-10 keV flux of SS433 as 8.7 x 10"" erg cm'^ s'^ which 
would correspond to a luminosity of 2.4 x 10^^ erg s~^ at the distance of M83. Thus, an 
SS433-like source would not necessarily be detected in X-rays with the current sensitivity. 

A second SNR with a decidedly bipolar morphology is M83-SNR-1-12 (BL37). This 
object also has spectoscopic data, and has strong [O III ] and [O II ] emission. 
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Fig. 1. — Continuum subtracted Ha images of the candidate SNRs found in the disk of M83. 
Each thumbnail is centered at the position of the candidate. The images have a square root 
stretch, and each covers 3x3 arcsec. on the sky. 
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3.2. Nuclear Candidates 

Identification of tlie SNR in tlie nuclear region is rather more difficult for several reasons. 
First, there is extensive diffuse interstellar emission which is also characterized by a high [S II 
]/Ha ratio. Second, the extreme reddening in some regions, coupled with the older underlying 
stellar population introduces an additional color term in the star subtraction process which 
leaves a residual in the [S II ] frame. Third, the starburst nature of the nucleus produces 
very strong and finely structured Ha emission, which makes the contrast between the H II 
regions and the SNR shells very low. 

As a consequence of these problems we used a different technique in the nucleus. Judging 
from the weakness of [O III ] A5007 emission in the whole of this field, we may safely assume 
(Kewley & Dopita 2002[ ) that the chemical abundances in the ISM are above, and possibly 



much above, the solar values. With such high chemical abundances, the temperatures of 
H II regions are quite low, possibly below 5000K, and at these temperatures collisional 
excitation of the [O II ]AA 3727,9 emission lines is effectively suppressed. In shock-excited 
gas, however, the cooling region in which we find the O II ion is characterized by a much 
higher temperature (~ 12000K), and the collisional excitation rate of the [O II ]AA3727, 3729 
emission lines is high. We therefore mainly rely on the strength of the O II emission to identify 
SNR candidates in the nuclear region. Because of the high and variable dust extinction in 
the nuclear region, the [O II ] line is heavily attenuated. As a consequence, all the nuclear 
candidates that lie on the far side of the dust mid-plane will have been missed. This implies 
that our completeness in this region cannot be appreciably better than 50%. The list of 
nuclear SNR candidates so identified is given in Table [3j These are listed by increasing 
RA and are named by galaxy, nature and field, and prefix (N) to identify them as nuclear 
candidates {viz. M83-SNR-1-Nxx), For reference, we also give the measured coordinates of 
the nuclear cluster in this Table. 

For those candidates which are not too affected by crowding with H II regions, we have 
measured their Ha fluxes using the digiphot aperture photometry routines in IRAF. Because 
(at the time of writing) the on-orbit photometric zero point of the F657N filter had not been 
determined, the transformation between electrons per second and flux was determined by 



comparison with the fluxes given in the Blair & Long (2004) catalog for the dozen candidates 



we have in common. In some cases, we find H II regions close to or within the boundaries 
of the SNR, and these objects were excluded from the fit. We estimate that the resultant 



fluxes are good to about ±25%. As in Blair & Long (2004) we make no attempt to correct 



the observed fluxes for the contamination by the [N II ]A6584 line. Accounting for this could 
reduce the measured flux by as much as 20 — 40%. However, an accurate determination of 
this correction will have to await spectrophotometric observations. 
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Images of the nuclear candidates are shown in [O II ]AA3727, 9 emission in Fig|2} Again, 
the images are shown on a square root stretch to bring out the fainter features. Each 
thumbnail covers 3x3 arcsec. on the sky. In the nuclear candidates the [O II ] emission 
is very strong compared to underlying H II regions, and this provides an excellent contrast 
between the SNR and the background. 



MSl-SMR-l-mi MSMNfi-t-NM MS3-5HfM4»4 M63 SNJM NOS M43-5MK-H«fi 




Mas-SNR-l-NOd Maa^^-l^NW Ma3-SI^R-1-N1& MS1-SNR<1-N11 Maa-Slffl-l-NIS MBa-SNR-l-MtS 




MS3'SHM4I14 MS^SHIM'tilS MS^IM^II? Mffl^SNR^-Nia MS^SNR-l-tllSi 



Fig. 2. — Continuum subtracted [O II ] images of the nuclear candidate SNRs. Each 
thumbnail is centered at the position of the candidate, is on a square root stretch, and 
covers 3x3 arcsec. on the sky. The small dark specks represent incompletely subtracted 
cosmic ray events, while the white specks are over-subtracted blue stars. 
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3.3. Extended O-strong nebulae 

As a result of the high chemical abundance in this field, [O III ] emission is generally 
weak. Therefore, every object with significant [O III ] emission is interesting in some way or 
another. The point sources with strong [O III ] are planetary nebulae (PNe), and we detect 



dozens of un-cataloged PNe in this field. Herrmann et al. (2008) have reported 241 PNe in 



M83, although the vast majority of these lie outside the WFC3 field. All five catalogued 
PNe that are in the WFC3 field have been detected as point sources in the [O III ] frame. 

Extended sources enhanced in [O III ] may be either Wolf-Rayet ring nebulae, or else 
young supernova candidates in which we directly see the ejecta of the supernova event as it 
passes through the reverse shock. It is also possible that shocks of more normal SNRs may 
appear strong in [O III ] emission for a brief time at the start of their radiative phase. 

Fairly strong extended [O III ] was detected from many of the SNRs listed by BL04 
because shock velocities above 100 km s~^ can still doubly ionize oxygen. However, we also 
might expect any young, core collapse SNRs (similar to Cas A in our Galaxy or 1E0102- 
7219 in the Small Magellanic Cloud) to be dominated by [O HI] emission lines. Indeed, the 
emission lines normally used to find ISM-dominated SNRs may not be present at all. 

For these reasons, we concentrated on inspecting any obviously extended sources of 
[O III ] emission in the WFC3 field that did not already correspond with known or newly 
identified SNR candidates. There are only a handful of reasonably bright [O III] sources, the 
five most interesting of which are listed in Table |4j This may be a heterogeneous sampling of 
objects, and their true identity will require follow-up spectroscopy to decipher, but several 
of the objects appear consistent with a possible 0-rich SNR interpretation. The continuum 
subtracted emission line images of these sources are shown in Figure |3| For each object 
there are four panels corresponding to their appearance in [O II ], [O III ], Ha and [S II ], 
respectively. We now discuss the properties of each of these sources in more detail. 

M83-SRC-1-01 is a small arc of unusually strong [O III ] emission that aligns with a faint 
diffuse patch of Ha. No enhancement in either [S II ] or [O II ] is seen and there are no 
strong stellar sources coincident with the region. The nature of this object remains unclear. 
M83-SRC-1-02 is a second case of a (very) faint, extended [O III ] nebula within a larger 
region of faint Ha emission. Its nature is unclear 

M83-SRC-1-03 is a clustering of 0-knots centrally located within an Ha and [S II ] 
shell which itself lies within a fainter, more extended nebulosity. Coincident stellar sources 
complicate the interpretation somewhat, but this object looks like a possible analogue of 
N132D in the Large Magellanic Cloud, which ~ 3200 years after the explosion has an outer 
~25 pc ISM-dominated shell surrounding inner 0-rich debris from the SN ( [Blair et al.||2000D . 
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If this analogy is correct, the outer shell in this case is roughly half as large as the one in 
N132D (0.36 X 0.48", or 8.0 x 10.6 pc), which could indicate a smaller precursor-blown 
cavity surrounding the SN. 

M83-SRC-1-04 is a curious bright [O II ] and [O III ] feature within a more extended 
region of Ha and [S II ] emission. It is by far the brightest object seen in the exposure 
using the WFC3 [O III ] filter. There is significant Ha emission coincident with the strong 
0-nebula, but [S II ] is weak. The 0-nebula appears as a partial shell open to the north, 
although a nearly coincident star or very tight cluster complicates the interpretation. In the 
V-band image the image is slightly fuzzy, which could be the result of the [O III ] emission 
being passed by the broad V-filter. This object could be either an 0-rich SNR or possibly a 
high-excitation WR ring nebula. 

M83-SRC-1-05 is an interesting, relatively bright extended object. In Ha, a half shell is 
clearly seen open to the west and brightest on the north and south sides, with a diameter 
of ~11 pc. Diffuse Ha emission fills the shell and spills out to the west and north. The full 
ring is more obvious in [O II ], and is present but very faint in [S II ] (certainly too faint to 
be indicative of shock activity). The northern and southern Ha bright spots are enhanced 
in [O III ] (especially in the north) although the rest of the shell can be seen. A bright 
continuum source is coincident not with the center of this shell but rather the northern edge, 
directly adjacent to the brightest [O HI] emission. This nebula is probably photoionized and 
perhaps a WR shell of some type, but it is clearly unusual due to the strength of the [O III 
] line compared with most emission regions in the field. 



Clearly follow-on spectroscopy is warranted for all of these interesting candidates. This 
is planned using the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope (Dopita 
eraL][2007|). 
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Fig. 3. — Emission line images of the objects with unusually strong and extended [O III ] 
A5007 emission. The key to the emission line corresponding to each image is given at the 
bottom right hand corner. The thumbnail images are each 3x3 arcsec on the sky. The 
lower three objects may either be 0-rich young SNR or else high excitation Wolf-Rayet ring 
nebulae. 



3.4. Other X-ray / Emission line Sources 



In order to see whether any compact sources had been missed, we inspected the region 
around each of the 23 Soria & Wu (2003) Chandra sources in the field that did not correspond 
to known or new SNRs found with the Ha - [S II ] - [O II ] technique. The most complicated 
regions of the nucleus were excluded from this search because of field crowding and source 
confusion. From this exercise, one source stood out as having a high probability of being an 
0-rich SNR. This is Soria & Wu (2003) source 70, located just to the east of the bright nuclear 
starburst region. Figure |4] shows a 3" region near this source (same panels as used above), 
where the white circle (1.5" diameter) is centered on the nominal Soria & Wu position; RA 
(J2000)= 13:37:01.28, Dec. (J2000)= -29:51:59.9. A relatively strong and compact nebula 
is found in both [O II ] and [O III ], slightly northeast of the nominal position. A slightly 
diffuse Ha nebula also visible at this position, but the object is very weak in [S II ]. Although 
the V-band is not shown, there is no appreciable continuum source at this position, and so 
a stellar residual or stellar confusion is not an issue. 

Table A.l of Soria & Wu (2003) indicates that source 70 is a soft X-ray source, with 
roughly 70 counts in the 0.3 - 1 keV band and 40 counts in the 1-2 keV band. This is 
consistent with expectations for a young, Cas A-like SNR. We identify the nebular object in 
Figure |4] as the likely optical counterpart of the X-ray source, and a bona fide 0-rich SNR in 
M83. The faint Ha emission might arise in the X-ray ionized nebula surrounding the blast 
wave. 
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Fig. 4. — This 4-panel figure shows a 3" region surrounding the position of Soria & Wu 
(2003) Chandra source 70. The white circle is 1.5" in diameter centered on the source 70 
position. Just northeast of the center is a nebular object with relatively strong [O II] (upper 
left) and [O III] (upper right) emission. A compact Ha source is also seen (lower left), but 
is much fainter in [S II] (lower right). This source is very likely a bona fide 0-rich SNR in 
M83. 
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3.5. The Search for the Remnant of SN1968L 

Of the six historical supernovae observed in M83 since 1923, only one, SN1968L, lies 
within our field of view. This event was observed in the western side of the nuclear region 
and was classified as a Type II, but its position is poorly known. Published positions are 
notoriously inaccurate (see discussion in BL04). We have therefore reviewed the literature 



and reconstructed the position based on the images provided by Wood & Andrews (1974). 
Scaling the WFC3 F555W image to approximate the appearance of Wood & Andrews Figure 
Ic, we can reconstruct the position of this object. We estimate this new position to be 
accurate to better than 1" and our derived coordinates should be good to the accuracy of 
the WCS solution used for the other SNR candidates. 

This region is marked on appropriately scaled WFC3 images of the nuclear region, see 
Figure |5} Inspection of these images reveals an emission nebula with a strong, probably 
unresolved core and a fainter, marginally resolved halo in [O III]. This has a coincident 
slightly extended [S II] nebula at about 0.2 arcsec. distant from our reconstructed position 
of the supernova event. No obvious Ha counterpart is present although the complexity of 
Ha emission and the faintness of the object against the strong Ha background could mask 
this. No obvious counterpart is present in the [O II ] waveband, but this could simply be 
due to heavy foreground dust attenuation. The reconstructed position for this object is RA 
13:37:00.433 and Dec. -29:51:59.65 (J2000). 

What might be expected for such a young SNR? If we assume the ejecta have been 
expanding at 10^ km s~^, it would still be less than 0.5 pc in diameter, which corresponds 
to <0.022". The unresolved nature of the core [O III ] source is thus consistent with a 
possible young supernova remnant identification. It is not entirely clear whether the [S II ] 
nebula is extended or whether the appearance is impacted by the continuum subtraction of 
nearby objects. We note, however, the similarity to Cassiopia A in our Galaxy where oxygen 
and sulfur dominate the appearance of the optical spectrum. The slightly extended halo in 
[S II ], and the corresponding [O III ] - emitting region may possibly arise from blast- wave 



photoionized pre-supernova stellar ejecta, similar to that seen in SN1987A (Fransson et al. 



1989 Crotts & Heathcote 1991). If placed at the distance of M83, this ring-like nebulosity 



would have a diameter of ~ 0.06", and be marginally resolved on these WFC3 images. 

In conclusion, we tentatively claim the object in the small circle in the Figure to be the 
recovered optical counterpart to SN1968L. A spectrum showing high [O III ] velocities would 
confirm this identification. 



Fig. 5. — A 3" field of WFC3 data centered on tlie expected position of SN1968L. North is 
at the top, east at the left. The panels include (upper left) V-band, (upper right) [O III ], 
(lower left) Ha, and (lower right) [S II ]. All images are shown with a square root scaling and 
all emission line images have been continuum-subtracted. The X marks the reconstructed 
SN1968L position as described in the text, and the large circle (1" radius) is the conservative 
upper hmit on its accuracy. The smaller (0.2" radius) circle marks the optical counterpart 
candidate for the SN, based on the appearance in the [O III ] image. This circle is reproduced 
in the other panels for reference. A faint stellar source hes just NE of the [O III ] source. 
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4. Interpretation of the Observations 

4.1. Cumulative Number : Diameter Relationships 

According to classical theory, a SNR evolves through a number of different phases in 
its passage to the point where it merges with the surrounding interstellar medium. These 



are summarized in Dopita & Sutherland (2003). In brief, the SNR first expands ballistically, 



R (X T until it has swept up sufficient mass to allow a reverse shock to propagate through 
the ejecta to heat it to very high temperature. The SNR then expands adiabatically in its 
Sedov phase R oc r^/^ until the outer blast wave shock becomes slow enough and old enough 
to start radiating away the energy stored in the hot gas. During the ballistic, reverse shock 
and Sedov phases the SNR is seen mostly at X-ray or radio frequencies, except for a short 
period when stellar ejecta passing through the reverse shock becomes prominent at optical 
wavelengths, especially in the lines of [O I ], [O II ], [O III ] and [Ne III ]. In the following, 
radiative, phase R oc r^/'' (McKee & Ostriker 1977). This is the phase when the SNR is 



most likely to be observed at optical wavelengths, and in this phase the abundances inferred 
from the radiative shock theory will reflect those of the surrounding ISM rather than the 
stellar ejecta. Finally, the energy contained in the hot plasma is radiated away, and the SNR 
enters its final momentum-conserving "snowplough" phase ( Oort|[l946 ) with R oc t^^^ before 
it turbulently merges with the general ISM. 

Since the M83 SNR are detected at optical wavelengths, they are either in the radiative 
phase, or else in the (ejecta dominated) reverse shock phase. In our [S II ] and [O II ] 
selected sample, none appear to be young ejecta-dominated SNR. Therefore to the degree 
that we can regard supernova as exploding at a constant stochastic rate, we would expect 
the cumulative number : diameter relation to follow A^(< D) oc D^/^. 

The observed distribution is shown in Figure [6j The SNR divide according to their 
environment - the nuclear group is consistent with being observed in free expansion, while the 
disk SNR generally follow the slope expected for SNR in their radiative phase. Superimposed 
on this slope are clear stochastic fluctuations due to the randomness of SNe. 

Why should the nuclear population be distinct from the SNR in the spiral arms? The 
most likely explanation for this is that there is a high porosity in the nuclear ISM. Given that 
it is a high-pressure region, the un-ionized part of the ISM is likely to be crushed to higher 
density, and the filling factor of the hot ISM is correspondingly larger. In sweeping up the 
hot tenuous ISM, the SNR blast wave would not be appreciably slowed until it encounters 
one of the dense H I clouds, driving a much slower radiative shock into it. At this point, 
it becomes very bright in optical emission lines. According to this scenario, the pre-shock 
density in the nuclear SNR should be much higher than in the disk SNR. We will show below 
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that this is indeed the case. 

For the nuclear SNR we can estimate both the ages and the supernova rate in the nucleus 
by assuming that in its free expansion phase the SNR expands at 5000-10000 km s^^. There 
are 11 SNR with a diameter less than 15 pc, with an age < 750 — 1500 yr. This implies 
a supernova rate of one every 70-140 yr. Incompleteness in our SNR catalog may almost 
double this estimated rate to one every ~ 50 yr. These estimates are consistent with the 
fact that one SN has been observed in this region since 1923. 



4.2. Derived SNR Parameters 



With the Ha photometry and the measured diameters, we can derive some interesting 
parameters for these SNR candidates. These parameters are summarized in Table |5j 

The total luminosity, L, can be derived directly from the Ha luminosity. For this, we 
use the distance given above and use the transformation L = 75LHa (for a shock velocity of 
200 km s~^. This transformation is derived from the shock models of Dopita &: Sutherland 



(1996 


) and 


Allen et al. 


(2008 



the total luminosity is weakly dependent on the shock velocity. Between 200 km s ^ and 500 
km s~^, the total luminosity relative to the Ha luminosity increases by a factor of two. 

For the M83 SNR we adopt a characteristic shock velocity, Vg, of 200 km s~^. This 
figure is chosen on the basis of that it is typical of the optical velocity dispersions observed 



' [e.g. 


ShuU & Hippelein 


Bilikova et al. 


2007 


). An 



can be derived from the fact that the [N II ] and Ha lines are clearly resolved in the low 



resolution spectroscopy of the M83 SNR by Blair & Long (2004). This limits the SNR 
expansion velocities to less than ~ 300 km s~^. An extreme upper limit on the shock 
velocity is placed by the condition that the cooling timescale for the shocked gas has to be 
shorter than the dynamical age of the SNR - otherwise the shocks would only be visible in 



X-rays. The cooling timescale for a radiative shock is given by Dopita & Sutherland (2003): 



Tcool ~ 200 



VlOO 



4.4 



ZriQ 



where Vioo is the shock velocity in units of 100 km s^^, Z is the chemical abundance in the 
pre-shock gas relative to solar, and uq is the pre-shock number density. In what follows, 
we will take Z = 1. In practice, Z is probably rather higher than solar. All of our SNR 
candidates are, by definition, in their radiative phase. Thus, if we can estimate no, we can 
then use equation [T] in conjunction with the dynamical age to estimate the maximum shock 
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Fig. 6. — The cumulative number: diameter relation for the M83 SNR. Here we distinguish 
the disk SNR population from the SNR discovered in the nuclear regions. The slope, 
m, of the relation for the disk SNR is consistent with these SNR being in their radiative 
phase (m = 7/2). The fluctuations around the slope reflect the stochastic timing of the 
supernova events. The nuclear SNR follow a slope which suggests that they are mostly 
in free expansion (m = 1) until the blast wave collides with a dense cloud in the ISM, 
propagating a radiative shock into it. 
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velocity consistent with the SNR being in its radiative phase. 

We may use the total luminosity to estimate the radiative lifetime of these SNR. The 
radiative lifetime is defined by: 

T"rad = sec, (2) 

where is the initial kinetic energy of the SNR ejecta. Here we take Eq ~ 10^^ ergs. We 
now use the area of the SNR shock defined by the mean optical diameter, as measured from 
the SNR image to estimate the mean surface brightness of the shock wave, S (erg cm~^ s~^). 

This radiative lifetime simply tells us the maximum time over which it can remain in the 
radiative phase before the stored explosion energy is lost as readiation. The radiative lifetime 
is correlated with the true age of the SNR, but not in any simple manner. Theoretically, the 
SNR has already passed through its free expansion and Sedov (adiabatic) phases before it 
reaches the radiative phase, so the actual age of the SNR depends on how long it has spent 
in each of these, and also on how long it has already spent in the radiative phase before 
being observed. 

Next, we use this surface brightness to estimate the pre-shock density, no, derived from 



radiative shock theory (Dopita & Sutherland 1996 Allen et al. 2008). For a fully radiative 



shock, the surface energy flux across the shock front is, 5* = pv 12. Transposing to extract 
the number density, and assuming a Helium abundance by number of 0.1, this equation can 
be written as: 

log (no/cm-3) = 2.88 + log 5 - 3 log vioo- (3) 

Note that the pre-shock density derived in this way is highly sensitive to the assumed shock 
velocity. Clearly, a measurement of the velocity dispersions of thes SNR would considerably 
assist in the accuracy of the derivation of the inferred parameters for individual SNRs. The 
pre-shock densities so derived lie in the range 0.6 < rio/cm"^ < 1700 (see Table These 



densities are consistent with the range of volume densities inferred by Heiner et al. (2008) at 
the surface of the giant molecular clouds in M83; 0.1 — 400 cm~^ for the region of M83 lying 
within i?25- Since the phase of the ISM which these densities trace is probably that of the 
warm neutral medium (WNM), the appropriate temperature is probably ~ 5 — 10 x 10'^ K, 
which implies that the range of ISM pressures in M83 is of the order 5 x 10'^ < P/k < 
1.5 X 10*^ cm-^K. 

Finally, as described above we use the density derived from equation |3] in conjunction 
with the radiative age given with equation |2] and the cooling timescale given in equation [T] 
to estimate the maximum shock velocity consistent with the requirement that the SNR is 
in the radiative phase. The values so determined are listed in Table [5j As expected, all of 
these are larger than the 200 km s~^ adopted as the likely shock velocity 
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4.3. The Radiative Age : Density Relationship 

Because the coohng timescale of the shocked gas scales inversely as the pre-shock den- 
sity, the distribution of the product rioTcooi should trace out the degree either to which local 
conditions influence the evolution (through variation in local density), or alternatively, the 
range in the explosion energy of the SNR candidates, Eo (which influences the radiative cool- 
ing timescale). The histogram of the product noTcooi is given in Figure [7j This distribution 
is clearly bimodal, suggesting that either the explosion energy is bimodal, or that there are 
strong and distinct environmental effects which form two populations of SNR. 

This bimodality can be investigated in a different way. We would also expect to see these 
SNR trace out a sequence with a slope of -1 on a logTcooi : log no diagram. This is plotted in 
Figure [sj Here we have divided the two populations at log [(no/cm~^)(rcooi/yr)] = 5.2 into a 
"high" and a "low" sample distinguished by yellow circles and the blue circles, respectively. 
These two distinct sequences are separated by a factor of four in density at constant age, or 
else a similar factor in terms of age at constant density. These two groups of SNR are not 
identical with the nuclear and the disk groups defined above, although most of the nuclear 
objects fall within the "high" sample. 

The likely measurement errors for the observable quantities which go into this rela- 
tionship are shown as error bars in Figure |8] It is clear that separation between these 
two sequences is too large to be explained by measurement errors. It therefore probably 
represents either that there are two precursor populations of Type II supernovae, giving 
explosion energies separated by a factor four, or that there is an environmental difference 
which provides a density contrast of about the same factor. 

In order to test the environmental hypothesis, we plot the SNR on the Ha image, 
distinguishing the "high" sequence points (yellow) from the "low" sequence ( marked in 
blue) - see Figure [oj For some of these SNR in locally crowded fields, it was not possible to 
measure the Ha fluxes. These are shown in green. From this figure it is evident that the 
high density sequence objects are nearly all confined either to the nuclear region, or to the 
spiral arm, while the low density sequence is more generally distributed in the disk. The 
most likely explanation for this extraordinary result is that the environment of the spiral 
arms and the nucleus provides higher density H I clouds for the SNR to interact with, most 
likely the result of a higher pressure in the interstellar medium in these environments. 

Since the sound speed in the X-ray emitting plasma is very high, of order 300-500 
km s~^, and the cooling timescale of this plasma is long, a local hydrostatic equilibrium can 
be achieved in this phase of the ISM. Therefore, in a medium with high porosity, the hot 
ISM can be thought of as supplying the substrate pressure confinement for the warm and 
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Fig. 7. — The histogram of the product noTcooi for our SNR candidates. The width of 
this should reflect both the range of local physical conditions in the ISM and the range 
of explosion energies of the objects. Remarkably, the histogram is clearly bipolar about 
log(norcooi) = 5.2, with the two peaks being separated by 0.5-0.6 in the logarithm of the 
product. 
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Fig. 8. — The relationship between the radiative age and the derived pre-shock density of 
the M83 SNRs. The error bars represent the displacement that would be produced by a 
±25% error in the measurement of the Ha flux, and by a ±10% error in measurement of the 
SNR radius. The theoretical slope of this relation is -1. It is clear that there are two distinct 
populations of SNR distinguished here by the different symbols. For a given dynamical age 
the upper (yellow points, high density) sequence has a pre-shock density four times higher 
than the lower (blue points, low density) sequence. 
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molecular phases of the ISM. The local pressure in the X-ray gas is in turn correlated with 
the surface brightness, provided that the scale height of the X-ray emitting plasma and the 
temperature of this plasma do not vary too much across the face of the galaxy. Thus, the 
surface brightness in X-rays should correlate with the local density, ^x-ray ^o- 

We have measured the surface brightness of the X-ray emission using the Chandra 



observations described by Soria & Wu (2003). These show diffuse emission concentrated 



toward the spiral arms and nucleus (in addition to the 127 point sources). The diffuse 
emission arises primarily from hot gas, as evidenced by the fact that spectra show strong 
features from lines expected in plasmas of 10^ - 10^ K. To determine the X-ray surface 
brightness at the position of each SNR, we created 0.4-2 keV images of the emission from 
M83, subtracting contributions from foreground diffuse X-ray emission and an extragalactic 
background. We measured the average surface brightness within 5" of each SNR, omitting 
portions of the circular regions affected by point sources. 

We compare the density derived from the optical data to the surface brightness (in 



arbitrary units) derived for the diffuse X-ray emission in Fig 10 The line drawn on this 
figure has the theoretical slope Sx-ray oc Kq. Scatter in this relationship is produced by 
variation in X-ray scale height. X-ray temperature, absorption of the soft X-rays as well as 
all the factors which enter into the determination of the pre-shock density in the SNRs. The 
absorption of soft X-ray emission in the ISM is more likely to be important in the nuclear 
regions, where the ISM column density may be > 10^^ cm~^. This affects the points with 
no > 500 cm~^. 



It is clear that there also is a bimodal behavior of the points plotted on Fig [TOj The 
"high" sequence points conform to the theoretically expected S'x-ray oc nl relationship, 
showing that the X-ray plasma does indeed supply the substrate pressure to the ISM in both 
the spiral arms and in the nucleus. The "low" sequence shows no such relationship. For 
these objects, the substrate pressure in the ISM must be provided by some other source, 
such as interstellar turbulence. The diffuse X-ray emission which is spatially coincident with 
these SNRs probably originates in the halo of the galaxy, with a much higher scale height 
than in the arm regions, and is not directly physically connected with conditions in the disk. 



The virtue of Fig 10 is that, unlike Fig|8| the bimodality cannot be ascribed to a bimodal 
behavior of an intrinsic property of the SNR population such as explosion energy Eq, but 
instead has to be related to the environment. 
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Fig. 9. — The positions of the M83 SNR are here plotted on the Ha image of M83. A 
number of the nuclear SNR are not shown because of crowding, and a lack of accurate 
photometry. Here we distinguish the objects by color. The yellow points represent the 
objects identified from figure [8] as belonging to the "high density" sequence, the blue points 
represent the SNR in the "low density" sequence, and the green points are for unclassified 
SNR. It is clear that the high density sequence objects are all confined either to the nuclear 
region, or to the spiral arm. 
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Fig. 10. — The correlation between the derived pre-shock density of the M83 SNRs and the 
local (5 " radius) diffuse X-ray surface brightness. Scatter in this relationship is produced 
by variation in X-ray scale height, X-ray temperature, absorption of the soft X-rays and 
extinction in the SNRs. The mean measurement error in the points is shown. If the X-ray 
plasma provides the substrate pressure in the ISM, and if scale height and the temperature 
of the X-ray plasma remained constant, then S'x-ray oc tiq. A line of this slope is plotted on 
the figure. It is clear that only the "high density" sequence (yellow points) conforms to such 
a relationship. Other factors are determining the pressure in the ISM for the SNR in the 
"low density" sequence (blue points). 
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4.4. Mass of Precursor Stars 



The mean mass and lower mass limit of the SNe precursor stars can be estimated from 
the number of observed SNR - as long as we know both the star formation rate and the mean 
age of the SNR at the completeness limit in diameter. The argument goes as follows. 

First, we can use the total flux in Ho; in the field to estimate the total star formation 
rate. The theory behind this is that the fiux in any hydrogen line is simply proportional 
to the number of photons produced by the hot stars, which is in turn proportional to the 
birthrate of massive stars. This relationship has been well calibrated at solar metallicity for 
the Ha line. In units of M0 yr~^, the estimated star formation rate is given by (Dopita & 



Ryder 


1994; 


Kennicutt 


1998 1 



[SFRuJMq yr-'] = (7.0 - 7.9) x lO'^^ [Lna/eTg.s~'] 



(4) 



For the full field of the continuum-subtracted WFC3 Ha image, we measure a total 
luminosity of Lhq- = 3.73 x lO^^erg s~^, which translates by equation |4] to a star formation 
rate of 2.76 Mq jt~^. This estimate does not take account of those H II regions which are 
obscured by dust, or of the unknown contribution of the [N II ] emission within the bandpass 
of the filter. 

The total number of SNR found in this image is 60. We may hope that any incom- 
pleteness in this number is matched by any incompleteness in our estimate of star formation 
rate, since both SNR and H II regions are observed at Ha and therefore subject to the same 
attenuation by dust. This assumption is certainly an over-simplification, but it cannot be 
tested in any simple way with our data. 

In order to compare the expected number of SNR with the observed number we need to 
know both the shape of the upper initial mass function (which determined how many stars 
are available to explode as SNe for a given star formation rate), and the mean (physical) 
age of the SNR at the point where the sample becomes incomplete. From the cumulative 
number : diameter relationship, it is apparent that the sample becomes seriously incomplete 
at a diameter of ~ 30pc at which point we have accumulated 46 SNR. We therefore attempt 
to estimate the age of the D < 25pc sample, which is the largest diameter at which we 
think the sample is complete. However, at this diameter we do not know the mean physical 
age, only the observed range of radiative lifetimes; 3.8 < log(ri.ad/yr) < 5.0. Assuming 
that the radiative lifetime is comparable with the expansion timescale at the time they are 
observed, we can estimate the true age of the SNR to be on average 2/7 of the observed 
radiative age. This (very crudely) puts the mean age of the SNR in the complete sample at 
log(r/yr) = 3.9 ± 0.5. 
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Once we have settled on the form of the IMF, we can demand that the number of SNRs 
we see up to the completeness limit matches the observed number of SNRs. The mean 
lifetime of the SNR population at the completeness limit provides an estimate of the lower 
mass limit of Type II SNe. The form of the IMF also provides the mean mass of SNe. We 
have adopted a Miller & Scalo (1979) IMF. This has 10.2% of the stars (by mass) with 
mass above 9 Mq and 3.6% above 20 Mq. On this basis, we estimate that the lower mass 
of Type II SNe in M83 is Mmm = ^(^-l Mq and that the mean Type II precursor mass is 

Mav = 24l^° Mq. 

Whilst these number are necessarily very approximate, we have demonstrated that the 
number of SNR detected in these images is consistent with expectations of stellar evolution 
theory and with the observed rate of star formation in the observed M83 field. 



5. Conclusions 

This WFC3 dataset has provided an rich source of physical insight into the nature of the 
supernova remnant population in M83, and their interactions with the interstellar medium. 
To summarize our results: 

• We have increased by a factor of five the number of SNR candidates in this field. 

• We have likely recovered the remnant of the historical supernova, SN1968L. 

• We discovered that the nuclear [R < 300 pc) SNR are physically distinct from those 
in the disk, being characterized by free expansion in a hot ISM with both high filling 
factor and pressure. We infer a supernova rate of around one every ~ 50 yr, consistent 
with the fact that one supernova event has been observed in this region since 1923. 

• The range in pre-shock density that we infer, 0.6 < no/cm~^ < 1700, is consistent with 



the Heiner et al. (2008) estimate of the densities at the surface of the giant molecular 
clouds of M83, obtained using GALEX far-ultraviolet (FUV) fluxes, VLA images of 
the 21-cm HI column densities, and estimates of the local dust abundances. 

There are two distinct populations of SNR divided at log [(r;,o/cm^^)(rcooi/yi')] = 5.2 
into "high density" and "low density" samples. The "high" sample is mostly confined 
to the nucleus and the spiral arms, and the "low" group is a field population. These 
populations apparently reflect the nuclear or arm and interarm ISM pressure regimes. 
For the nuclear region and the arm region, the substrate pressure is provided by the 
diffuse X-ray plasma. For the inter-arm region, the substrate pressure is not provided 



by the X-ray gas, and is more likely provided by turbulent motion in the molecular 
and warm phases of the ISM. 
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Table 1. Log of WFC3 Observations (Prop ID# 11360). 



Filter 


Exposure (s) 


Date (2009) 


F336W 


600 


August 26 


F373N 


3 X 800 


August 20 


F438W 


3 X 640 


August 26 


F438W 


10 


August 26 


F487N 


3 X 900 


August 25 


F502N 


3 X 828 


August 26 


F555W 


3 X 401 


August 26 


F555W 


10 


August 20 


F657N 


4 X 371 


August 25 


F673N 


2 X 600 


August 20 


F673N 


650 


August 20 


F814W 


3 X 401 


August 26 


F814W 


10 


August 26 



Table 2. Catalogue of Outer {R > 300pc) SNR Candidates 



Name RA (J2000) Dec. (J2000) D(arcsec)i D(pc) F(Ha) ^ F(Hq;)b&l^ B&L# 



M83-SNR-1-01 


13 


36 


58.712 


-29 


51 


00.52 


1.14 


25.2 


l.lE-14 








X41, St 


M83-SNR-1-02 


13 


36 


58.885 


-29 


52 


26.02 


0.92 


20.4 












M83-SNR-1-03 


13 


36 


59.166 


-29 


51 


47.96 


0.56 


12.4 


2.4E-15 








X66 


M83-SNR-1-04 


13 


36 


59.479 


-29 


52 


03.66 


0.68 


15.0 


1.7E-14 


1.64E- 


■14 


BL31 


X45 


M83-SNR-1-05 


13 


36 


59.725 


-29 


50 


32.88 


1.18 


26.1 


1.7E-15 










M83-SNR-1-06 


13 


37 


00.068 


-29 


52 


08.43 


0.52x0.86 


15.3 


9.0E-15 










M83-SNR-1-07 


13 


37 


00.343 


-29 


51 


20.63 


0.96x1.32 


25.2 


7.0E-15 


4.20E- 


•15 


BL33 


Strong 


M83-SNR-1-08 


13 


37 


00.672 


-29 


52 


21.75 


0.96x1.44 


26.5 












M83-SNR-1-09 


13 


37 


01.022 


-29 


50 


56.31 


0.88 


19.5 


3.6E-15 










M83-SNR-1-10 


13 


37 


01.076 


-29 


51 


41.63 


0.56 


12.4 


4.9E-15 










M83-SNR-1-11 


13 


37 


01.579 


-29 


52 


04.93 


1.12 


24.8 












M83-SNR-1-12 


13 


37 


01.729 


-29 


51 


13.47 


1.40 


31.0 


2.1E-14 


2.05E- 


•14 


BL37 


X77,biF 


M83-SNR-1-13 


13 


37 


02.002 


-29 


52 


50.03 


2.2 


48.7 


2.1E-15 


7.32E- 


•15 


BL38 


small H 


M83-SNR-1-14 


13 


37 


02.115 


-29 


51 


58.75 


0.64x2.54 






1.16E- 


■14 


BL39 


Jet- like 


M83-SNR-1-15 


13 


37 


02.443 


-29 


51 


26.07 


0.56 


12.4 


1.9E-14 


2.28E- 


■14 


BL41 


X81, im 


M83-SNR-1-16 


13 


37 


02.631 


-29 


51 


09.23 


1.08 


23.9 










small H 


M83-SNR-1-17 


13 


37 


04.877 


-29 


52 


18.59 


1.28 


28.3 


1.8E-15 










M83-SNR-1-18 


13 


37 


05.807 


-29 


52 


46.01 


1.68 


37.2 


1.7E-15 








Very fa: 


M83-SNR-1-19 


13 


37 


06.435 


-29 


50 


25.19 


1.58 


34.9 


1.4E-14 


1.44E- 


■14 


BL47 




M83-SNR-1-20 


13 


37 


06.986 


-29 


51 


09.45 


0.62 


13.7 


1.6E-15 










M83-SNR-1-21 


13 


37 


07.686 


-29 


51 


09.93 


0.38 


8.4 












M83-SNR-1-22 


13 


37 


07.471 


-29 


51 


33.23 


1.32x1.58 


32.1 


2.3E-14 


1.46E- 


•14 


BL53 


X105 


M83-SNR-1-23 


13 


37 


07.595 


-29 


52 


19.63 


1.28 


28.3 


4.1E-15 


5.80E- 


•15 


BL55 




M83-SNR-1-24 


13 


37 


07.852 


-29 


52 


41.93 


1.26 


27.9 


1.4E-15 








R47? 


M83-SNR-1-25 


13 


37 


08.110 


-29 


52 


21.31 


0.56x1.08 


18.1 


3.3E-15 










M83-SNR-1-26 


13 


37 


08.324 


-29 


50 


56.43 


1.40 


31.0 












M83-SNR-1-27 


13 


37 


08.490 


-29 


52 


02.35 


2.30 


50.9 










Bright ] 


M83-SNR-1-28 


13 


37 


08.559 


-29 


51 


34.87 


0.74 


16.4 


1.4E-14 


1.19E- 


■14 


BL58 


X109 


M83-SNR-1-29 


13 


37 


08.667 


-29 


52 


42.69 


1.26 


27.9 


6.1E-15 










M83-SNR-1-30 


13 


37 


08.656 


-29 


51 


53.57 


0.72 


15.9 


5.9E-15 










M83-SNR-1-31 


13 


37 


08.696 


-29 


51 


33.29 


1.84 


40.7 


1.3E-14 










M83-SNR-1-32 


13 


37 


08.753 


-29 


51 


37.41 


1.26 


27.9 


8.6E-15 


1.18E- 


■14 


BL59 




M83-SNR-1-33 


13 


37 


09.044 


-29 


51 


33.43 


0.64 


14.2 


3.0E-15 








R52? 
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Tablc 2 — Continued 



Name 


RA (J2000) 


Dec. (J2000) 


D(arcsec)^ 


D(pc) 


F(Ha) 2 


F(Ha)B&L^ 


B&L# 


M83-SNR-1-34 


13 


37 


09.227 


-29 


51 


34.06 


0.60 


13.3 


2.5E-15 






M83-SNR-1-35 


13 


37 


09.318 


-29 


50 


58.50 


1.24 


27.4 


3.9E-14 






M83-SNR-1-36 


13 


37 


10.117 


-29 


51 


28.22 


0.48 


10.6 


1.6E-14 






M83-SNR-1-37 


13 


37 


10.329 


-29 


51 


28.74 


1.14 


25.2 








M83-SNR-1-38 


13 


37 


10.381 


-29 


51 


34.18 


1.18 


26.1 


2.7E-15 






M83-SNR-1-39 


13 


37 


10.797 


-29 


51 


44.85 


1.82 


40.3 


2.8E-15 


5.08E-15 


BL60 


M83-SNR-1-40 


13 


37 


10.833 


-29 


52 


44.52 


0.68 


15.0 


1.4E-15 







Notes^ 
R52? 



^Typical measurement error ±0.08 arcsec. 

^Units: erg cm"^ s~^, typical measurement error ±15%. 

3Prom Table 5 of Blair & Long (2004). 

^X: X-ray source number from Soria &; Wu (2003); R: radio source number from Maddox et al. 
(2006) using VLA data. 



Table 3. Catalogue of Nuclear {R < 300 pc) SNR Candidates 



Name 


RA (J2000) 


Dec. (J2000) 


D(arcsec)^ 


D(pc) 


F(Ha) 2 


Notes' 


M83 Nucleus 


13 37 00.871 


-29 51 55.97 










M83-SNR-N-01 


13 37 00.039 


-29 52 02.08 


0.20 


4.4 






M83-SNR-N-02 


13 37 00.204 


-29 51 58.47 


0.66x0.90 


17.3 






M83-SNR-N-03 


13 37 00.214 


-29 52 06.27 


0.68 


15.0 


4.9E-14 


X53 


M83-SNR-N-04 


13 37 00.335 


-29 52 05.52 


0.32 


7.1 


2.0E-14 


X54 


M83-SNR-N-05 


13 37 00.373 


-29 51 58.88 


0.24x0.14 


4.2 




Near SN1968L 


M83-SNR-N-06 


13 37 00.385 


-29 52 01.72 


0.84 


18.6 






M83-SNR-N-07 


13 37 00.404 


-29 52 06.33 


0.52 


11.5 






M83-SNR-N-08 


13 37 00.545 


-29 51 58.88 


0.56 


12.4 


4.1E-13 




M83-SNR-N-09 


13 37 00.581 


-29 52 09.24 


1.84 


40.7 






M83-SNR-N-10 


13 37 00.609 


-29 51 59.65 


0.76 


16.8 


3.9E-13 


X59? 


M83-SNR-N-11 


13 37 00.661 


-29 51 57.10 


0.35x0.2 


6.1 






M83-SNR-N-12 


13 37 00.701 


-29 51 59.89 


0.86 


19.0 




X59? 


M83-SNR-N-13 


13 37 00.757 


-29 52 06.49 


0.44 


9.7 




X62 


M83-SNR-N-14 


13 37 00.929 


-29 51 54.25 


0.88 


19.5 




X63 


M83-SNR-N-15 


13 37 00.932 


-29 51 55.93 


0.18 


4.0 






M83-SNR-N-16 


13 37 01.107 


-29 51 52.13 


0.26x0.68 


10.4 




X67 


M83-SNR-N-17 


13 37 01.172 


-29 51 57.19 


0.52 


11.5 






M83-SNR-N-18 


13 37 01.317 


-29 51 57.17 


1.20 


26.5 






M83-SNR-N-19 


13 37 01.593 


-29 52 02.29 


0.96 


21.2 




X74 



Typical measurement error ±0.08 arcsec. 
^Units: erg cm~^ s~^, typical measurement error ±15%. 

^X: source number from Soria &: Wu (2003); source crowding and coordinate uncertainties make 
some of these identifications tentative. 



Table 4. Extended [O III] Nebulae in M83 WFC3 Field 1 



Name 



RA (J2000) Dec. (J2000) D(arcsec) D(pc) 



Notes^ 



M83 Nucleus 
M83-SRC-1-01 
M83-SRC-1-02 
M83-SRC-1-03 
M83-SRC-1-04 
M83-SRC-1-05 



13 37 00.871 
13 37 04.120 
13 37 06.976 
13 37 07.132 
13 37 09.058 
13 37 09.761 



-29 51 55.97 
-29 51 03.85 
-29 50 57.14 
-29 51 14.27 
-29 52 09.92 
-29 52 43.93 



0.45x0.20 10x4 [O III] arc in diffuse Ha. 

0.33 7.3 0-knots within faint Ha; possible 0-ric 

0.3x0.5 6.7x11.1 [O HI] knots within Ha ring; 0-rich SN 

0.36 8.0 Young 0-rich SNR in H H region? 

0.50 11.1 Strong [O HI] and Ha ring; possible W 



^R: source number from Maddox et al. (2006) (VLA data). 



Table 5. Derived Parameters of M83 SNR. 



iNctine 


1 J 


r yncx ) 


lr>rr T 2 
log Li 


iOg Trad 


log D 


log no 






M83-SNR-1-01 


25.2 


l.lE-14 


39.31 


4.19 


-0.95 


1.04 


L 


460 


M83-SNR-1-02 


20.4 
















M83-SNR-1-03 


12.4 


2.4E-15 


38.65 


4.85 


-0.99 


0.99 


H 


640 


M83-SNR-1-04 


15.0 


1.7E-14 


39.50 


4.00 


-0.31 


1.68 


H 


580 


M83-SNR-1-05 


26.1 


1.7E-15 


38.50 


5.00 


-1.79 


0.20 


L 


450 


M83-SNR-1-06 


15.3 


9.0E-15 


39.22 


4.28 


-0.60 


1.38 


H 


580 


M83-SNR-1-07 


25.2 


7.0E-15 


39.11 


4.39 


-1.14 


0.84 


L 


460 


M83-SNR-1-08 


26.5 


2.2E-14 


39.61 


3.89 


-0.69 


1.30 


L 


450 


M83-SNR-1-09 


19.5 


3.6E-15 


38.83 


4.68 


-1.21 


0.78 


H 


520 


M83-SNR-1-10 


12.4 


4.9E-15 


38.96 


4.54 


-0.68 


1.30 


H 


640 


M83-SNR-1-11 


24.8 
















M83-SNR-1-12 


31.0 


2.1E-14 


39.59 


3.91 


-0.84 


1.14 


L 


420 


M83-SNR-1-13 


48.7 


2.1E-15 


38.59 


4.91 


-2.24 


-0.25 


L 


340 


M83-SNR-1-14 


















M83-SNR-1-15 


12.4 


1.9E-14 


39.55 


3.95 


-0.09 


1.89 


L 


640 


M83-SNR-1-16 


23.9 


1.7E-15 


38.50 


5.00 


-1.71 


0.27 


L 


470 


M83-SNR-1-17 


28.3 


1.8E-15 


38.52 


4.98 


-1.83 


0.15 


L 


440 


M83-SNR-1-18 


37.2 


1.7E-15 


38.50 


5.00 


-2.09 


-0.11 


L 


390 


M83-SNR-1-19 


34.9 


1.4E-14 


39.42 


4.09 


-1.12 


0.86 


L 


400 


M83-SNR-1-20 


13.7 


1.6E-15 


38.47 


5.03 


-1.25 


0.73 


H 


610 


M83-SNR-1-21 


7.5 


1.2E-15 


38.33 


5.17 


-0.87 


1.10 


H 


800 


M83-SNR-1-22 


32.1 


2.3E-14 


39.63 


3.87 


-0.84 


1.15 


L 


410 


M8S-SNR-1-2S 


28 3 


4.1E-15 


38 88 


4.62 


-1 47 


n 51 


L 


440 


M83-SNR-1-24 


27.9 


1.4E-15 


38.42 


5.09 


-1.93 


0.05 


L 


440 


M83-SNR-1-25 


18.1 


3.3E-15 


38.79 


4.71 


-1.18 


0.80 


H 


540 


M83-SNR-1-26 


31.0 
















M83-SNR-1-27 


50.9 
















M83-SNR-1-28 


16.4 


1.4E-14 


39.42 


4.09 


-0.47 


1.52 


H 


560 


M83-SNR-1-29 


27.9 


6.1E-15 


39.05 


4.45 


-1.29 


0.69 


L 


440 


M83-SNR-1-30 


15.9 


5.9E-15 


39.04 


4.46 


-0.82 


1.17 


H 


570 


M83-SNR-1-31 


40.7 


1.3E-14 


39.38 


4.12 


-1.29 


0.69 


H 


470 


M83-SNR-1-32 


27.9 


8.6E-15 


39.20 


4.30 


-1.14 


0.84 


L 


440 


M83-SNR-1-33 


14.2 


3.0E-15 


38.75 


4.75 


-1.01 


0.97 


L 


600 
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Table 5 — Continued 



Name 


D 


F(Ha) ^ 


1 T 2 

log L ^ 


log Trad 


log S ^ 


1 5 

log no ° 


Class ^ 


^ \ 7 

t>s(maxj 


M83-SNE-1-34 


13.3 


2.5E-15 


38.67 


4.83 


-1.03 


0.95 


H 


620 




27 4 


3.9E-14 


39 86 


3 64 


-0 47 


1.52 


L 


440 






1.6E-14 


39 47 


4 03 


-0 03 


1 95 


H 


680 


M83-SNE-1-37 


25.2 


1.6E-14 


39.47 


4.03 


-0.78 


1.20 


L 


460 


M83-SNR-1-38 


26.1 


2.7E-15 


38.70 


4.80 


-1.59 


0.40 


L 


450 


M83-SNR-1-39 


40.3 


2.8E-15 


38.72 


4.78 


-1.95 


0.04 


L 


370 


M83-SNR-1-40 


15.0 


1.4E-15 


38.42 


5.09 


-1.39 


0.59 


L 


580 


M83-SNR-1-N02 


17.3 


l.OE-14 


39.27 


4.23 


-0.66 


1.32 


H 


550 


M83-SNR-1-N03 


15.0 


4.9E-14 


39.96 


3.54 


0.15 


2.14 


H 


580 


M83-SNR-1-N04 


7.1 


2.0E-14 


39.56 


3.94 


0.41 


2.39 


H 


820 


M83-SNR-1-N07 


11.5 
















M83-SNR-1-N08 


12.4 


4.1E-13 


40.88 


2.62 


1.24 


3.23 


H 


640 


M83-SNR-1-N10 


16.8 


3.9E-13 


40.85 


2.65 


0.95 


2.93 


H 


560 


M83-SNR-1-N16 


10.4 
















M83-SNR-1-N18 


26.5 
















Meas. Error ^ 


±1.8 


±10% 


±0.15 


±0.15 


±0.20 


±0.20 




±25% 



1 —2 —1 

erg cm s . 
■^Total luminosity: erg s~^. 
^Radiative age: yr. 

^Mean surface energy flux of radiative shocks: erg cm~^ s~^. 
^Inferred pre-shock density: cm~^. 
^Density class: high-H, low-L. 
"^Maximum shock velocity: km s^^. 

'^Errors given here are typical, and based upon measurement errors only. Errors in the assump- 
tions used may produce much larger systematic errors. These can be estimated from the formulae 
given in section [4?2| 



